
THE AEROSPACE 

CORPORATION 

El Segundo, California 



AJASfl-CA- R(e/3C 

AEROSPACE REPORT NO. 
ATR-93(7251 )-1 


Analytic Model of Aurorally Coupled 
Magnetospheric and Ionospheric 
Electrostatic Potentials 


15 January 1994 NASA-CR-196130 

19940030852 


Prepared by 
J. M. CORNWALL 

Space and Environment Technology Center 
Technology Operations 


Prepared for 
NASA 

Goddard Space Flight Center 
Greenbelt,MD 20771 


Grant No. NAGW-2126 


Engineering and Technology Group 


PUBLIC RELEASE IS AUTHORIZED 



LANGLEY RESEARCH CENTER 
LIBRARY NASA 
HAMPTON, VIRGINIA 










TECHNOLOGY OPERATIONS 


The Aerospace Corporation functions as an "architect-engineer" for national security 
programs, specializing in advanced military space systems. The Corporation’s Technology 
Operations supports the effective and timely development and operation of national security 
systems through scientific research and the application of advanced technology. Vital to the 
success of the Corporation is the technical staffs wide-ranging expertise and its ability to stay 
abreast of new technological developments and program support issues associated with rapidly 
evolving space systems. Contributing capabilities are provided by these individual Technology 
Centers: 

Electronics Technology Center: Microelectronics, solid-state device physics, VLSI 
reliability, compound semiconductors, radiation hardening, data storage 
technologies, infrared detector devices and testing; electro-optics, quantum 
electronics, solid-state lasers, optical propagation and communications; cw and 
pulsed chemical laser development, optical resonators, beam control, atmospheric 
propagation, and laser effects and countermeasures; atomic frequency standards, 
applied laser spectroscopy, laser chemistry, laser optoelectronics, phase conjugation 
and coherent imaging, solar cell physics, battery electrochemistry, battery testing and 
evaluation. 

Mechanics and Materials Technology Center: Evaluation and characterization of 
new materials: metals, alloys, ceramics, polymers and their composites, and new 
forms of carbon; development and analysis of thin films and deposition techniques; 
nondestructive evaluation, component failure analysis and reliability; fracture 
mechanics and stress corrosion; development and evaluation of hardened 
components; analysis and evaluation of materials at cryogenic and elevated 
temperatures; launch vehicle and reentry fluid mechanics, heat transfer and flight 
dynamics; chemical and electric propulsion; spacecraft structural mechanics, 
spacecraft survivability and vulnerability assessment; contamination, thermal and 
structural control; high temperature thermomechanics, gas kinetics and radiation; 
lubrication and surface phenomena. 

Space and Environment Technology Center: Magnetospheric, auroral and cosmic 
ray physics, wave-particle interactions, magnetospheric plasma waves; atmospheric 
and ionospheric physics, density and composition of the upper atmosphere, remote 
sensing using atmospheric radiation; solar physics, infrared astronomy, infrared 
signature analysis; effects of solar activity, magnetic storms and nuclear explosions 
on the earth’s atmosphere, ionosphere and magnetosphere; effects of electromagnetic 
and particulate radiations on space systems; space instrumentation; propellant 
chemistry, chemical dynamics, environmental chemistry, trace detection; 
atmospheric chemical reactions, atmospheric optics, light scattering, state-specific 
chemical reactions and radiative signatures of missile plumes, and sensor out-of- 
field-of-view rejection. 



DISPLAY 94N35358/2 
94N35358*# ISSUE 11 CATEGORY 46 

RPT# : NASA-CR-196130 NAS 1.26:196130 ATR-93‘(7251)-1 CNT# : NAGW-2126 

94/01/15 12 PAGES UNCLASSIFIED DOCUMENT 

UTTL: Analytic model of aurorally coupled magnetospheric and ionospheric 
electrostatic potentials 

AUTH: A/CORNWALL, J. M. 

CORP: Aerospace Corp., El Segundo, CA. CSS: (Engineering and Technology Group. 

) 

SAP: Avail: CASI HC A03/MF A01 

CIO: UNITED STATES 

MAJS: /*AURORAL ARCS/*AURORAL ZONES/*AURORAS/*ELECTROSTATICS/*IONIZATION/*POLAR 
REGIONS 

MINS: / ANALYSIS (MATHEMATICS)/ DISSIPATION/ MATHEMATICAL MODELS/ NONLINEARITY/ 
PHENOMENOLOGY 

ABA: Author (revised) 

ABS: This paper describes modest but significant improvements on earlier 

studies of electrostatic potential structure in the auroral region using 
the adiabatic auroral arc model. This model has crucial nonlinearities 
(connected, for example, with aurorally produced ionization) which have 
hampered analysis; earlier work has either been linear, which I will show 
is a poor approximation or, if nonlinear, either numerical or too 
specialized to study parametric dependencies. With certain simplifying 
ENTER: MORE 



DISPLAY 94N35358/2 

assumptions I find new analytic nonlinear solutions fully exhibiting the 
parametric dependence of potentials on magnetospheric (e.g.. cross-tail 
potential) and ionospheric (e.g., recombination rate) parameters. No 
purely phenomenological parameters are introduced. The results are in 
reasonable agreement with observed average auroral potential drops, 
inverted-V scale sizes, and dissipation rates. The dissipation rate is 
quite comparable to tail energization and transport rates and should have 
a major effect on tail and magnetospheric dynamics. This paper gives 
various relations between the cross-tail potential and auroral parameters 
(e.g., total parallel currents and potential drops) which can be studied 
with existing data sets. 

ENTER: 



Aerospace Report No. 
ATR-93(7251)-1 


ANALYTIC MODEL OF AURORALLY COUPLED 
MAGNETOSPHERIC AND IONOSPHERIC 
ELECTROSTATIC POTENTIALS 


Prepared by 
J. M. Cornwall 

Space and Environment Technology Center 
Technology Operations 


15 January 1994 


Engineering and Technology Group 


THE AEROSPACE CORPORATION 
El Segundo, CA 90245-4691 


Prepared for 
NASA 

Goddard Space Flight Center 
Greenbelt, MD 20771 


Grant No. NAGW-2126 


PUBLIC RELEASE IS AUTHORIZED 



Aerospace Report No. 
ATR-93(7251)-1 


ANALYTIC MODEL OF AURORALLY COUPLED 
MAGNETOSPHERIC AND IONOSPHERIC 
ELECTROSTATIC POTENTIALS 


Prepared 



Approved 



J. B. Blake, Director 
Space Sciences Department 



A. B. Christensen, Principal Director 
Space and Environment Technology 
Center 


NOTE 


The material reproduced in this report originally appeared in Journal of Geophysical Research. 
The ATR is published to document the work for the corporate record. 



JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 98, NO. A9, PAGES 15,373-15,379, SEPTEMBER 1, 1993 


Analytic Model of Aurorally Coupled Magnetospheric and Ionospheric 

Electrostatic Potentials 

John M. Cornwall 1 


Space and Environment Technology Center, The Aerospace Corporation, El Segundo, California 

This paper describes modest but significant improvements on earlier studies of electrostatic 
potential structure in the auroral region, using the adiabatic auroral arc model. This model has crucial 
nonlinearities (connected, for example, with aurorally produced ionization) which have hampered 
analysis; earlier work has either been linear, which I will show is a poor approximation or, if nonlinear, 
either numerical or too specialized to study parametric dependencies. With certain simplifying 
assumptions I find new analytic nonlinear solutions fully exhibiting the parametric dependence of 
potentials on magnetospheric (e.g., cross-tail potential) and ionospheric (e.g., recombination rate) 
parameters. No purely phenomenological parameters are introduced. The results are in reasonable 
agreement with observed average auroral potential drops, inverted-V scale sizes, and dissipation rates. 
The dissipation rate is quite comparable to tail energization and transport rates and should have a 
major effect on tail and magnetospheric dynamics. This paper gives various relations between the 
cross-tail potential and auroral parameters (e.g., total parallel currents and potential drops) which can 
be studied with existing data sets. 


1. Introduction 

In this paper I look again at some issues raised over the 
years concerning the adiabatic auroral arc model [Chiu and 
Schulz, 1978; Lyons, 1980, 1981; Chiu and Cornwall , 1980; 
Chiu et al. f 1981]. The essence of this mode is a linear 
relation between the auroral parallel current J|| and the 
potential drop iff = tpj — tp e between the ionospheric 
electrostatic potential ip{ and the equatorial potential *p e 
along an auroral field line (see equation (2) in section 2). In 
principle, this model can account for all the gross features 
(e.g., potential drops, auroral size scales, and dissipation 
rates) of magnetosphere-ionosphere coupling as governed by 
auroral phenomena, with no adjustable phenomenological 
parameters or fitting factors, as was pointed out by the 
author some years ago [ Cornwall , 1983]. However, because 
the auroral model is nonlinear beyond the current-potential 
relation, it has not yet been possible to give, even in an 
idealized model, a precise analytic picture of how potential 
drops and so forth depend on physically determinate param¬ 
eters such as the polar cap potential drop, ip pc and the 
ionospheric recombination rate constant, a. The studies 
referred to above were either linear or treated nonlinearities 
numerically for the most part (although Cornwall [1983] did 
give some nonlinear scaling laws which will be refined here). 

In later studies [ Cornwall , 1988, 1990], certain exact 
solutions were found to the nonlinear model equations, but 
these were too specialized to allow a study of parameter 
dependence; typically, some combination of independent 
parameters had to be fixed to allow for a solution. The 
present investigation is carried out in a similar spirit but with 
different assumptions and results, which allow most of the 
interesting physical parameters to be varied freely. In par¬ 
ticular, I model the equatorial potential tp e as Vasyliunas 
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[1970] did, which simplifies the nonlinear structure at the 
price of introducing unphysical discontinuities on the bound¬ 
ary between open and closed field lines. These discontinui¬ 
ties are of little importance to my major results; their 
primary effect is to make unreliable the detailed shape of 
if/, - if/ e very near the central auroral field line. 

For present purposes one essential nonlinearity of the 
adiabatic auroral arc model is associated with a density- 
dependent Pedersen conductivity. (This and other nonlinear¬ 
ities contribute on smaller size scales to auroral structure 
and instabilities [ Cornwall , 1990; Keskinen et aL , 1992], 
which will not concern us here.) It is a strong nonlinearity in 
the sense that the ionospheric plasma density on the central 
auroral line can be 10 or more times greater than the 
unperturbed density. It can, of course, be dealt with numer¬ 
ically, as the early studies referred to previously in this 
section did, but that is not my purpose here. We will find that 
it significantly affects the gross size of auroras, increasing 
their width by a factor of several over those of a linearized 
analysis [e.g., Chiu et aL, 1981; Cornwall, 1983] to a value of 
several hundred kilometers. A second nonlinearity arises in 
the relation between the auroral potential drop ip and the 
height-integrated ionospheric plasma density, N, which re¬ 
sults in nonlinear relations between \p and tp pc .and between 
the dissipated auroral power P and ip pc (the appropriate 
scaling laws were given earlier by Cornwall [1983]). 

Our major results on the relations between the auroral 
size, ip, P , and ip pc are given in equations (43H49); with no 
adjustable parameters, they give values within a factor of 2 
or better of observed values for typical values of *p pc . 

These results have interest not only in themselves but also 
for the future studies they suggest. Section 1 concludes with 
some remarks in this direction, setting the present work in a 
larger context involving tail and magnetospheric convection 
dynamics. 

The standard picture of magnetospheric convection is 
largely based on the work of Vasyliunas [1970] (see also 
Fejer [1964]; Schield et al. [1969] for important precursors). 
In this picture, magnetic field lines are equipotentials, and 
the two-dimensional electrostatic potential, mapped to the 
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ionosphere, is governed by ionospheric current conserva¬ 
tion. The effective height-integrated conductivities receive 
contributions both from the ionosphere and from the ring 
current, and there are parallel currents driven by the con¬ 
vection discontinuity at the boundary between open and 
closed field lines. This picture has been modeled with 
considerable success in computer simulations [e.g., Wolf 
1970; Wolfet a /., 1991]. 

Of course, field lines in the auroral zone are not equipo- 
tentials, as many studies have shown f Evans , 1974; Mizera 
and Fennell , 1977; Reijf et al., 1988; Lindqvist and Mark - 
land, 1990]. The convection modelers have not yet added the 
physics discussed in this paper to their codes, in part 
because it adds substantial complexity and, perhaps, in part 
because it seems at first glance to have impact localized to 
the auroral zone. 

I suggest that the analytic treatment of this paper could be 
used, at least at first, as a substitute for adding complex 
codes to the convection models. My results could be used as 
phenomenological input to existing model codes. I further 
suggest that some treatment of auroral dissipation will be 
essential, if the modelers are ever to hope to capture the 
main features of tail and magneto spheric dynamics, even in 
regions remote from the auroral zone. The reason is that 
auroral-ionospheric Joule dissipation is [ Cornwall , 1983] 
comparable to energization and transport rates in the tail and 
ring current and so plays an important role in global energy 
balance. (A similar point was made by W. Liu, as cited by 
Huang et al. [1989].) 

Consider the following simple argument. The energy per 
unit ionospheric area stored in a dipolar flux tube is 
~L 4 R E P , where R E is the radius of the Earth and P-4 x 
10“ 9 ergs cm' 3 is a typical plasma sheet pressure, while the 
Joule dissipation rate 7L p E\ is of the order of 10 ergs cm' 2 
s or more during an aurora. So the time scale r on which an 
aurora could drain this flux tube of energy is 


l 4 r e p 


T » 


1 P E 


2 

X 


~ 10 2 - 10 


( 1 ) 


not large compared to tail transport and energization times. 
(Note that the tail energy is stored mostly in ions, while the 
aurora transfers energy from the magnetosphere to the 
ionosphere mostly via energetic electrons.) It has been long 
known [Erickson and Wolf, 1980; Schindler and Birn , 1982] 
that there are serious difficulties in creating a steady state 
loss-free model of tail convection. This suggested the possi¬ 
bility of nonsteady convection as discussed by the above 
authors or lossy convection (as indicated, for example, by an 
effective adiabatic index of y < 1 [Spence et al. , 1987]). An 
important contributor to loss is cross-tail drift [Kivelson and 
Spence , 1988; Spence and Kivelson , 1990], and equation (1) 
suggests that ionospheric Joule dissipation may be just as 
important. 

I am currently in the process of constructing semianalytic 
models of the effect of auroral dissipation on tail dynamics, 
which may be of some use. Ultimately, I hope that the 
computer modelers will take over and vastly improve my 
simple efforts. 


2. Finding the Electrostatic Potential 

It is well known [ Vasyliunas , 1970] that the usual sort of 
two-dimensional electrostatic potential in and near the polar 


cap calls for nonvanishing 7|, driven by a nonvanishing 
divergence of horizontal ionospheric currents. Later, Corn - 
wall [1983] identified these parallel currents with auroral 
currents according to the standard prescription of the adia¬ 
batic auroral model: 

Ji=-G(^-(2) 

where *p i € is the ionospheric or equatorial potential along a 
given field line and Q is a parameter of the order of 
ne 2 (m e V e ) 1 with n as the plasma sheet electron density 
and V e as a characteristic plasma sheet electron velocity. 
This led to the prediction [ Cornwall , 1983] that 

*~<4-6)—^ (3) 

k e 

where 

A = (1 p IQ) m ~ 100 km (4) 

is the characteristic length scale of the adiabatic auroral arc 
model and i J/ pc is the polar cap potential drop. (Actually, (2) 
holds only when ijs pc is large enough to give ^ > 1 kV; for 
smaller </ahe necessary can be furnished by other means). 

The earlier work was largely linear, and even the linear 
analysis was not carried to completion. In this section I will 
give the full linear analysis and also an exact nonlinear 
solution to the model equations. In some respects the linear 
analysis is not a good approximation to the exact solution (it 
has the wrong spatial scale), but it can be plausibly fudged to 
give global results (e.g., the central potential drop) quite 
similar to the exact analysis. Moreover, it illustrates how I 
turn the equatorial potential ijj e , which appears as a source 
term in the nonlinear equations, into a boundary condition. 

The two fundamental equations of the adiabatic auroral 
arc model are 

Vj_ • (2 • = -7| = Q<p (5) 

dN r , 

— + V E -V ± N = -Q+-a(N 2 -N 2 0 ) (6) 

o i e 

where N is the height-integrated ionospheric plasma density, 
2 is the height-integrated ionospheric conductivity tensor, T 
is the number of electron-ion pairs produced per incident 
auroral electron, a is a height-integrated dissociative recom¬ 
bination coefficient, aN 2 summarizes nonauroral sources of 
ionospheric ionization, and V E is the electric drift velocity 
such that 

(7) 

The return current region lying just outside the auroral 
region is governed by different physics, which I will not 
attempt to model in detail here (for an earlier attempt, see 
Chiu et al. [1981]). Presumably, i/f — goes through zero at 
the edges of the aurora and then becomes slightly negative in 
the return current region, but with \if/ — \p e \ much smaller 
than in the auroral region. In the simplified analysis I will 
give, tf/ — {f/ e approaches zero only at asymptotically large 
distances from the zero, and so the return current is also far 
away. I emphasize that this is a mere artifact, which unless 
it were mistakenly taken literally, has no effect on our 
conclusions. In reality, the return current region lies just 
adjacent to the auroral region, only a few auroral scale 
lengths. A, away. 
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Equations (5-7) contain a vast number of effects which are 
readily studied with a computer but whose analytic treat¬ 
ment is either difficult or impossible. These include various 
instabilities, as well as nonlinearities, the primary effect of 
which is to increase the density N and conductivity 2 in the 
auroral region. This conductivity enhancement has a number 
of important and well-known consequences for convection, 
which will be studied elsewhere in connection with auroral- 
magnetospheric coupling. For now, I want to find an analyt¬ 
ically treatable treatment of this conductivity nonlinearity in 
the spirit of earlier [ Cornwall , 1988, 1990] analytic solutions 
to similar auroral problems. These earlier solutions are not 
useful in the present context, because they do not allow for 
a priori specification of the input driving potential tp e . 
Fortunately, i p e is not completely arbitrary, and progress 
can be made by modeling i p e as Vasyliunas [1970] and 
Cornwall [1983] did; tp e consists of two separate solutions to 
Laplace’s equation (reflecting charge neutrality), joined on a 
closed contour representing the auroral zone and continuous 
across this contour. However, V ip e is discontinuous across 
the contour. This discontinuity leads to J| proportional to a 
Dirac delta function. The input potential \p e -then only 
appears as a boundary condition on a homogeneous nonlin¬ 
ear equation, some special solutions to which can be found 
analytically. 

Now for the details. I begin by simplifying the ion con¬ 
ductivity equation (6), dropping the terms on the left-hand 
side. These terms are primarily important [i Cornwall , 1990; 
Keskinen et aL , 1992] for auroral instabilities on time scales 
of ~100 s, but on longer time scales it is a reasonable 
approximation to assume rough equilibrium between recom¬ 
bination and auroral precipitation. Then (6) yields 

( 8 ) 

and (5) is rewritten in terms of ip: 

V x • (2 • - Q<p = —V ± • (2 • V ± *,). (9) 

Under the assumption that the various parameters (a, T, 
N 0 ) which appear in (8) are constants, the Hall term drops 
out of the left-hand side of (9), since then both 2 and \p are 
functions of the single function N. Thus (9) becomes 

V x • (2 P V X( /,) - Qip = -V x * (2 • (10) 

Now we need to specify the source i p e and the right-hand 
side of (10). In the equipotential case (\p = 0) the right-hand 
side vanishes almost everywhere except on the convection 
boundary between open and closed lines, where it is a delta 
function. I give a simple example of such a source term, 
based on constant 2 and ip f as two different solutions of 
Laplace’s equation, joined at the convection boundary. Use 
spherical coordinates with the magnetic dipole axis as the z 
axis and the azimuthal angle <p defined to be tH2 at dawn and 
— W2 at dusk; 8 is the magnetic colatitude. Let 

j (tan j 8 ) sin <p 

'f'e = ~'i' P c -:- e<e pc (li) 

2 tan i 8 pc 

1 sin tp pc sin <f> 

'i'c^Z'f'pc -r-r- it ~ 8 pc > 8> 8 pc . (12) 

2 y sin 8 y 


(For the southern polar cap, replace 8 with n - 8 and 
tan i 6 pc with cot x - 6 pc in (11). Here \p pc is the dawn-to-dusk 
potentii across the polar cap, and the convection boundary 
is at 8 = 8 pc (or 0 = tt — 8 pc in the south). Both (11) and 
(12) satisfy Laplace’s equation in the angular variables. The 
form (11) represents an essentially constant dawn-dusk 
electric field over the polar cap, as one sees from the 
stereographic projection of the unit sphere from the south 
pole to the jc - y plane tangent to the north pole, with 

tan | 8 —» p; x = p cos <f> 

(13) 

<t> -* <P\ y = p sin <j>. 

Then (11), in the projected variables, is ip y, yielding a 
constant field. 

I will now make the idealization in (10) of saving only the 
most singular term, which is a delta function. It is this which 
allows further progress to be made, since the nonlinear terms 
on the left-hand side of (10) depend on ip e only through 
boundary conditions at 8 - 6 pc , it — 8 pc . The result is that 
(10) becomes, in the northern hemisphere. 


Vi* (2 p V.i» - Q*P = 


*PpcZp sin <p 
~2RlsinT‘ 


• (1 + cos d)8(e - e pc ) y (14) 

where, on the right-hand side, everything depending on 8 
(including 2 P ) is to be evaluated at 6 = $ pc . I will not 
explicitly write out the contributions from the southern 
hemisphere, which are easily supplied by symmetry. 

Note that in the idealization of saving only the most 
singular source term, the specific forms (11) and (12) for \p e 
and the assumption of constant 2 when 8 ¥=■ Q pc are 
irrelevant; all that matters is the coefficient of the delta 
function in (14). 

To find a solution to (14) two solutions to the homoge¬ 
neous version of this equation must be found, matched in 
value at 6 = 6 pci and the discontinuity in 6 derivatives 
adjusted to match the delta function. The two solutions are 
chosen so that each decays exponentially with angular 
distance from $ pc . Far from 8 = 8 pc the conductivity 2 P 
approaches, it is assumed, a constant value, at least in the 
sense that 2 P varies more slowly (owing, for example, to 
day-night effects) than it does due to precipitation in the 
auroral zone. In that case, (14) is a linear equation, straight¬ 
forwardly solved. Later in this section the linear version of 
(14) is solved, including the delta function source; this both 
yields the necessary linear solutions far from the auroral 
zone and illustrates the matching procedure used to accom¬ 
modate the delta function. 


2.1. Linear Case 

The equation to be solved for constant 2 p is 

i a / a«A l aV , 

—-sin a —) + ~r~ 2 —:-A 2 

sm 8 38 y 88) sm 2 0 sin 8 

sin 4> 

= K 8(8 - 0 pe ), (15) 

sm 8 


where 
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K - \ <l>pc( 1 + COS e pc ) 

A = R e A = R^Q!2 P ) m . 


( 16 ) (In the limit N a —> 0 the scaling law t p c ~ ip p< ? was already 
given by Cornwall [1983].) We will compare this equation 

( 17 ) later to the nonlinear solution developed next. 


Since A, the width of the (idealized) auroral zone, is about 

100 km, A ^> 1. We seek a solution of the form 2.2. Nonlinear Case 


*P = K sin <t>G(e, 0 pc ), (18) 

where G( 0 , $') is a Green’s function for Legendre’s equa¬ 
tion with m = 1, /(/ + 1 ) = -A 2 . That is, 

/= 4 ± I (1 “ 4A 2 ) 1/2 — “ 7 ± /A. (19) 

In the limit of large A the exact solution for G can be formally 
expressed as a linear combination of the conical functions 
^-i/ 2 +/a( cos 0)» 2-!-i/2+ja( cos 0) [ e -8*» Abramowitz and 
Stegun , 1964]. However, it is difficult to find useful forms for 
these rather obscure special functions, and we choose in¬ 
stead to give an approximate solution in terms of the 
well-known functions K x> I x (Hankel and Bessel functions 
of imaginary argument). For small 0 the homogeneous solu¬ 
tions to (15) are &i(A0) and /j(A0), valid when A0 3 1. 

Since A0 3 1, these do not cover a useful range in 0. A 

better approach is to seek homogeneous solutions of the 
form, for example, F(6)K X [A/?(0)], determining F and R so 
that dangerous terms in the differential equation proportional 
to A or A 2 are exactly cancelled. A straightforward compu¬ 
tation shows that this requires 


The nonlinearities are important only in the immediate 
vicinity of the auroral zone, which has a small width com¬ 
pared to R e . Therefore I will use the flat earth approxima¬ 
tion in the equations, with the x coordinate perpendicular to 
the auroral zone and y along it. The inhomogeneous delta 
function such as in (14) need not be written explicitly, since 
its only role is to furnish a boundary condition at the auroral 
zone. Equation ( 8 ) is used to eliminate ip in favor of N , 
assuming that N c and a/TQ in this equation are constants. 
With the observation that 2 P is linear in N y so that 32 p ldN 
is also constant, the homogeneous version of ( 10 ) is 

Vj. • (NV ± N 2 ) - Q(d2p/dN) -\N 2 -Nl)= 0 (24) 

or 

3 

Vi « 3 ~—1 (n 2 -1) = 0 (25) 

2Ag 

, Ipo „ N 

A 2 = -jr, Ipo = Zp(N = N„), . (26) 

** O 


( 0 V' 2 

R=e,F=[ — } (20) 

ysin 0 ] 

With this choice the fractional derivation of F^|(A/?) from 
solving (15) is 0 2 /36 for small 0 and 0(A~ 2 ) for large 0. 

It is easy to find the appropriate linear combinations to 
match the delta function in (15), with the final result that the 
solution to the linear problem of (15) is 



K sin <p 
2 


00 


pc 


sin 0 sin 6 pc 


1/2 

[*i(A 0 )/j(a e pc )H(e 


- V) + /i(a0)a: 1 (a e pc me pc - 0)], (21) 

where H{x) is the unit step function ( H{x ) = 1, x > 0; = 
0, x > 0). This potential falls off exponentially (roughly as 
exp [-A|0 - 0 pc |]) on both sides of 6 pc . Note that A 6 pc 
1, so that in the vicinity of the auroral zone one should use 
large-argument asymptotic expansions to evaluate the 
Bessel functions. When this is done, one finds for the 
potential drop at 0 = 0 pct sin <)> = -1: 


«/r c = <A(0 - 0 pc ) ~ 


A ip 


pc 


4R e sin 0 


— AE l . 


( 22 ) 


pc 


In extending this linear solution to a nonlinear solution the 
point is that A depends on and thus on the height- 
integrated plasma density N, which by ( 8 ) depends on ip. A 
simple-minded approach to this nonlinearity is to use ( 8 ) to 
express N (and 2 p ) in terms of ip, and then (22) becomes an 
algebraic equation for tp c : 

1/2 / r«^"i .t. \ 1/4 




sin <pip 


pc 


1 02 , 


4R e sin 0 pc \Q dN 



TQiP^ 


ae 


Nonlinear equations similar to (24) have been studied before 
analytically [ Cornwall , 1988, 1990] and numerically [Keski- 
nen et al ., 1992] with full two-dimensional dependence. 
However, the known two-dimensional analytic solutions are 
too restricted for the present purpose, and my present 
strategy is to save in (24) and (25) only gradients in the x 
direction (across the aurora, roughly north-south). Thus I am 
not modeling auroral blobs and transient (i.e., unstable) local 
structures. In any event, in the absence of structural mag- 
netospheric forcing the gradient structure which is more or 
less persistent is the x direction gradient which is saved. 
There are also examples [ Cornwall , 1988] of Kelvin- 
Helmholtz stable fully two-dimensional structures in which 
the x and y direction length scales are essentially the same 
and roughly equal to the scale length I will find. 

So I will ignore the y variation and replace (25) with 

f'-|a; 2 (f 2/3 -i) = o, (27) 

where F — n 2 and prime indicate dldx. This equation can be 
reduced to quadratures; I quote the solution in terms of the 
original variable n: 

n c dnn 2 / 2\ 3 10 

, (^T) [r + 3j + 27 

This integral supplies the boundary conditions that at x = 0, 
n — n c = N c /N 0 , where N c is the central ionospheric 
density in the auroral zone; typically, n c » 1 . Another 
constant of integration has been supplied so that as \x\ -+ °°, 
n —> 1 at an exponential rate: 

n -> 1 + C x exp (-xlA 0 ) + 0[exp (-2x/A 0 )]. 



- 1/2 


5 1/2 A, 


( 28 ) 


(23) 


X —► <» 


(29) 
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We will soon see that the constant Cj in (29) is determined by 
n Ci which is in turn determined by matching a discontinuity 
in n’(x = 0) to as in the linear case. 

The integral in (28) can be reduced to a combination of 
integrals of rational functions plus an elliptic integral of the 
third kind [e.g., Whittaker and Watson , 1952] and so can be 
said to yield an exact solution in terms of known functions, 
at least in principle. But this form is of little use, because it 
gives x in terms of n while we want n in terms of x. 
Nevertheless, for the dedicated reader who wishes to pursue 
it, I quote the change of variables which yields a standard 
elliptic integral form: 

n + f+± 10 1/3 (1 - 3 I/2 ) 

u - -. (30) 

« + | + i 10 ,/3 (l + 3 1/2 ) 

To determine the constant C \, one may use (29) to replace 
xlA 0 in (28) to find 



and the upper sign in (35) is for x < 0 and the lower sign for 
x > 0. This form of the solution is useful out to a|x| = 1 
-« C “ I/2 , where n - 0(1). 

For large \x\ a similar analysis leads to 

n = 1 + C,e ±x,A '- (7/6)[C,exp (±x/A 0 )] 2 + (101/48) 

• (C|e =x,A ') 3 + • • • (37) 

with the same sign convention as in (36). Here C\ is the same 
constant as in (32). Note that it is qualitatively (but not 
necessarily quantitatively) accurate to estimate Cj by find¬ 
ing a point Xj where the slopes and values of the large \x\ and 
small |x| forms of n in (35M37) match. Just saving the first 
nontrivial terms yields 

*1 = (20) m A 0 (nl n - Cj' 2 ) = (20) m (A c - C\ n A„) 



(38) 

C\ = (C 2 - 1) exp (X|/A 0 ) 

(39) 

1 


C 2 = — (11 + 21 1/2 ) = 1.56. 

10 

(40) 


+ In (n - 1). (31) 


By adding and subtracting an integral, which can be explic¬ 
itly done and which removes the singularities in the inte¬ 
grand of (31) at both endpoints, one can find 


where 




exp [2(/3 - 1) + y] 


(32) 



(33) 





= y(n c = eo) + - (5 ,/2 )«; I/2 + 0(«;') 


-i 


(34) 


Tc(«c = ®) = -2.50 

where the last form of (34) is useful for large n c . 

Rather than continue with analyzing the integral (28), one 
finds a more useful form of the solution for large and small 
x\ by direct solution of the differential equation. For small 
x| this amounts to an expansion in powers of /i c -2 <$: 1 and 
leads to 


n = /i c (l ± ax) 


1 +-*; 2 (1 ± ax)~ A 


50 

567 


n c 4 (1 ± ax) 8 + • • • 


(35) 


It follows that for large n c (either from the above equations 
or from the more accurate (32)), the transition from the 
nonlinear regime (40) to the linear regime (saving only the 
first term of (37)) sets in at a value of x of the order of 5A C , 
where A c is the auroral scaling length CL P /Q) 1/2 evaluated at 
the center of the auroral zone. In effect, x j is the distance 
from the center of the aurora to the unperturbed ionosphere, 
and is several times larger than one would have supposed 
from a linear analysis. 

The scaling xj — (20) 1/2 A,. is one of several key results 
following from the nonlinear analysis. The second result is a 
formula for the central density N c or, equivalently, for if/ c = 
ifi(x = 0) using (8). It is derived as follows: First, the 
inhomogeneous term for the nonlinear equation (27) is sup¬ 
plied from the fundamental equation (14), using x = Re(Q pc 
- 6) and the relation between ^ and N in (8). The result for 
the duskside aurora where sin <f> = -1 is 




V - 


2A: 


(* 2 -l)=- 


3H2 


2aeNll po 


'I'pcZpcil + cos 0 PC ) 


pc> 


2 R e sin 6 pc 


\S(x) 


One inserts the solution (35) to find 


„ 5 125 

nan c \\ + -n c - — n c + 


—A5(x). 



(41) 


(42) 


It turns out that n c 1, so only the first term of this 
expansion is saved. Then, using the definition of a in (36) and 
rearranging, one finds 


5 1/2 r / 3S p \ 1/2 » pc (l + cos V) 
8 ae dNJ R E sin 6 pc 


(43) 


where 


a = [2A 0 (5n c ) m ]~ l 


Putting in nominal ionospheric numbers (T = 100, a = 10 14 
cm 2 sec -1 , Q = 0.1 cm -1 sec -1 , and d2 p /dN = 3 cm 3 

(36) 

sec -1 , one finds 
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2/3 


pc 


60 kV, 


(44) 


(I use 60 kV as a reference polar cap potential, because a 
similar value is associated with the threshold for substorm 
occurrence [see Weimer et al. , 1992; Ahn et al., 1992].) 

The third key result expresses the central potential drop 
if/ e in terms of the polar cap potential drop, Q, and iono¬ 
spheric parameters and is just a combination of (8) and (43): 


ae 



(N 2 C - N 2 0 ) = 


aeN 2 

~TQ~ 


2/3 



Ippcd + COS 0 DC ) 


PC' 


16 \Qae 
Using the same constants as before 

tpc 


Re sin 6 pc 


4/3 


(45) 


4/3 


if/ c = 3.4 kV 


160 kV 


(46) 


This result can be compared to the N 0 = 0 version of (23), an 
equation based on a simple-minded fudging of the linear 
theory. One can easily check that (23) in the N 0 = 0 limit 
gives (45) except that the constant factor 5 2/3 =*0.183 in (45) 
is replaced by 4” 4/3 0.157. The agreement is remarkably 

close. 

A final key result is the auroral dissipation of power, 
integrated over the auroral zone. I will define a theoretically 
useful, if not immediately observationally relevant, excess 
dissipation power P as that due to the field E ± - -V ± if/ 
associated with the difference (if/ = if/ { - if/ e ) between the 
ionospheric potential with and without auroras. So 


P = 2 R e sin 6 pc J dx d<f> 2 p (V ± iIf) 2 (47) 

where the factor of 2 counts both polar caps and the integrals 
over dx (distance across the aurora) and d<f> extend over the 
auroral zone. To be definite, I will use for the x dependence 
of N and if/ the first term in (35), valid for rt c » 1, integrate 
in x over —a “ 1 < x < a “ 1 , assume that N and if/ depend on 
sin <f> as their dependence on if/ pc would suggest, that is, if/ - 
|sin <f>| 4/3 and N ~ |sin ^l 273 , and integrate <f> over the region 
sin <f> ^ 0. The result is 


P =(—) r( * )r( ^ } ** c(1+cos 8pc)1 
\ 64 / T(j) \ ae /\ aAr / r e sin 2 0 pc 

(48) 

P =10” ergs sec- • (49) 

This is substantial dissipation, quite enough to influence 
substorm and tail processes as discussed in connection with 
(1) in section 1. 

I am unaware of any studies in the literature which deal 
directly with the main results in (44), (45), and (46). But it is 
clear that the predicted numbers are, in an average sense, in 
reasonable agreement with auroral observations [e.g., 
Evans, 1974; Re iff, 1988; Lindqvist and Marklund , 1990; 
Chiu et al., 1982]. 



Fig. 1. Sealed density n - N/N 0 versus scaled length ax (for a, 

see (36)), for n c = 10. 

3. Concluding Remarks 

My main results are (36), giving the inverted-V scale 
length a~ l as V20A C , where A c is the linear scale length 
(see (4)), based on the central Pedersen conductivity, and 
(43M49), giving N c , if/ c , and the dissipated power P. The 
scale length is rather bigger than one would have estimated 
from the linear theory in section 2, and so the linear profiles 
are not very good approximations. However, the ad hoc 
equation (23), which grafts certain elements of the nonlinear 
theory on to the linear analysis, gives a value for if/ c in good 
agreement with the full nonlinear analysis. This analysis 
gives values for N c , if/ c , and P (with no arbitrary parame¬ 
ters), which are in decent agreement with observed values, 
and yields scaling laws for these quantities in their depen¬ 
dence on ip c which can be experimentally studied. 

One cannot expect fully realistic spatial profiles of ^and N 
near the center of the auroral zone, because my fundamental 
hypothesis for the magnetospheric potential if/ e (equations 
(11) and (12)) yields cusps in if/ and N at 6 = 6 pc or x = 0. 
These are illustrated in Figures 1 and 2, in which N and if/ are 
plotted for a typical value of n c = 10. In principle, N and if/ 
should approach x = 0 with zero slope, since these are 
symmetric around x = 0. There is no reason to distrust these 
profiles when the density and potential have dropped by, 
say, a factor of 2, and the dissipated power P is an integral 
over the profiles which is insensitive to the cusp behavior. 
The cusps will, of course, be removed by a smoother choice 
of if/ e , but precisely how to make this choice is not presently 
clear. 

Having analytic and fully parametrized expressions for, 
for example, the dissipated power P (sec (48)) will be very 
useful for theoretical investigations of tail transport and 
dynamics. As mentioned earlier, this power is comparable to 
that lost in cross-tail drift, and can have a fundamental 
impact on the picture of disturbed tail phenomena. Investi¬ 
gations to be reported later are now underway on this 
subject. 

It is, of course, important to know whether the predictions 
of this paper concerning the relation between if/ pc and 
various auroral phenomena are borne out by data. The 
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Fig. 2. Ionosphere-equatorial potential difference scaled by the 
central value (for which, see (45) and (46)), for n c - 10. 

material necessary for such a study, in the form of convec¬ 
tion field data (yielding 4/ pc ), field-aligned current data, and 
field-aligned potential drops, all exists and needs to be 
correlated. The intervals for which all these data are simul¬ 
taneously available is by no means negligible. A study of this 
sort would be of great interest quite aside from the present 
theoretical considerations. 
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